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ABSTRACT 

We report on Keck and Hubble Space Telescope observations of the binary millisecond 
pulsar PSR B1855+09. We detect its white-dwarf companion and measure mpsssw — 
25.90 ± 0.12 and mF8i4w = 24.19 ± 0.11 (Vega system). From the reddening-corrected 
color, (mFsssw — "^F8i4w)o = 1-06 ± 0.21, we infer a temperature T e s = 4800 ± 800 K. 
The white-dwarf mass is known accurately from measurements of the Shapiro delay of 
the pulsar signal, Mq = 0.258to'xm M Q . Hence, given a cooling model, one can use the 
measured temperature to determine the cooling age. The main uncertainty in the cooling 
models for such low-mass white dwarfs is the amount of residual nuclear burning, which 
is set by the thickness of the hydrogen layer surrounding the helium core. From the 
properties of similar systems, it has been inferred that helium white dwarfs form with 
thick hydrogen layers, with mass ;> 3 x 10~ 3 M Q , which leads to significant additional 
heating. This is consistent with expectations from simple evolutionary models of the 
preceding binary evolution. For PSR B1855+09, though, such models lead to a cooling 
age of ~ 10 Gyr, which is twice the spin-down age of the pulsar. It could be that the 
spin-down age were incorrect, which would call the standard vacuum dipole braking 
model into question. For two other pulsar companions, however, ages well over 10 Gyr 
are inferred, indicating that the problem may lie with the cooling models. There is no 
age discrepancy for models in which the white dwarfs are formed with thinner hydrogen 
layers (^3 x 1O~ 4 M ). 

Subject headings: binaries: general — pulsars: individual (PSR B1855+09) — stars: 
evolution — white dwarfs 
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1. Introduction 

About one twentieth of the known radio pulsars 
reside in binary systems. For most of these, the 
companions have estimated masses between 0.1 and 
0.4 M© and are thought to be low-mass, helium white 
dwarfs. Generally, the properties of the pulsars in 
these systems differ markedly from those of typical 
isolated pulsars, showing more rapid spin periods and 
smaller inferred magnetic fiel ds (for reviews, see, e.g ., 
Phinncy fc Kulkarni 1994| ; |Van den Heuvel 19951) . 



Presumably, this is because of a phase of mass and 



an gular momentum tranafcr which occurred when the 



with mass ;> 10 -3 Mq, residual nuclear burning in the 
p-p cycle will be a significant source of heat during 
the further evolution, and, for a given age, the white 
dwarf will be hotter than one whose hydrogen layer 
is thinner. 

Evidence in favor of thick hydrogen layers and 
significant nuclear burning comes from the millisec- 
ond pulsar binary system PSR J1012+5307. For this 
binary, ipsR — 7 Gyr (Lorimer et al. 1995), while 
from the surface temperature and gravity, as deter- 
mined using optical spectroscopy (T ff ~ 8500 K, 
logy ~ 6.7; Van Kerkwijk, Bergeron, & Kulkarni 



pro genitor of the current white dwarf aacended on the 

giant branch and overfilled its Roche lobe. At the ces- 
sation of mass transfer, the neutron star turned on as 
a millisecond radio pulsar and the companion was left 
as a helium white dwarf. 



1996; Callanan, Garnavich, & Koester 199? ), one in 



An interesting property of these binaries is that 



fers a much shorter cooling age, tvvD — 0-5 Gyr, if 
one uses models in which the white dwarf is assumed 
to have formed with a relatively thin hydrogen layer 
(^3 x 10 - 4 M Q ), with, in consequence, little nuclear 
burning (Lorimer et al. 1995 ; Barna, Marks, fc Con 



the y contain two independent clocks that started run- 



ning more or less simultaneously: the spin-down age 
of the pulsar and the cooling age of the white dwarf 
(Kulkarni 1986). Assuming that the pulsar spins 
down due to a torque N oc v n , the pulsar age is given 

by 



non Smith 1996; Althaus & Benvenuto 1997; Hansen 



fc Phinncy 1998| ). It has been argued that this im- 



ipsR = 



(n - 1)P 



1 - 



(1) 



where P = 1/v is the current spin period, P is its rate 
of change, Pq is the period when the pulsar began 
spinning down following cessation of mass transfer, 
and n — vv/v 2 is the "braking index," equal to 3 
under the assu mption of magnetic dipolc radiation 
(for a review, see Lyne fc Smith 1998 ). For n — 3 and 
Po <C P, ipsR — t c = P/2P, where r c is the pulsar 
"characteristic age." 

The cooling age of the white dwarf, the second 
clock, can be determined from the stellar tempera- 
ture and mass using a cooling model. The cooling 
properties of helium white dwarfs have been modeled 
extensively, particularly since the optical identifica- 
tion of a number of white-dwarf companions of mil- 
lisecond pulsars. A central issue that has arisen is 
how much hydrogen remains when the white dwarf is 
formed. Most likely, the amount of hydrogen left will 
be anti-correlated with core mass, since for larger core 
masses the pressure at the core-envelope interface re- 
quired for CNO cycle shell burning can be maintained 
down to lower envelope masses. If a white dwarf is 
left with a sufficiently thick remaining hydrogen layer, 



plies that the pulsar has not yet spun down much, 
i .e., P p ~ P. Alberts, Sa vonije , & Van den Heuvel 
( 1996 ) and Driebe et al. ( 1998| ), however, find from 
simple evolutionary models in which the mass loss 
during the red giant phase is simulated, that the he- 
lium white dwarf companion should have had a much 
thicker hydrogen layer at formation, ~5 x 10~ 3 M©. 
With this thicker layer, they obtain a cooling age that 
is consistent with the pulsar spin-down age. 

The results are less clear for other systems with op- 
tically identified companions, since these are too faint 
for optical spectroscopy and their masses can only be 
estimated from the mass function (using a guess for 
the pulsar mass and statistical arguments for the or- 
bital inclination). There is one system, however, for 
which the mass of the white dwarf is known precisely. 
This system is PSR B1855+09, composed of a 5.4 ms 
radio pulsar and a low-mass companion in a 12.3d 



circular orbit (Segelstein et al. 1986). The pulsar is a 



very stable rotator and it has been possible to measure 
the general relativistic Shapiro delay of the pulsar sig- 



nal near superior conjunction ( Ryba fc Taylor 1991 



Kaspi, Taylor, fc Ryba 1994 ). From the measure 



mcnts, one infers that the binary is nearly edge on 



(sini = 0.9992 



+0.0004 
-0.0007 



) and that the companion has 



f0 028 M Q . In addition, from the 
1.1 ± 0.3 mas, the distance 



O.258_ ool6 iw 



7T 



mass Mq 
measured parallax 
is 0. 91^^20 kpc. The characteristic age r c of the pul- 
sar is 5 Gyr. Optical observations of the field have so 
far failed to detect the companion, but set stringent 
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limits, showing that it must be an old, cold white 
dwarf ( Callanan ct al. 1989| ; Kulkarni, Djorgovski, & 
Kl dmola 1991^ 

Here, we report on Keck and Hubble Space Tele- 
scope observations of the PSR B1855+09 held, in 
which the counterpart is detected. We use these to 
determine the effective temperature and discuss the 
implications for cooling of helium white dwarfs and 
braking of millisecond pulsars. 

2. Observations 

The PSR B1855+09 held was observed on 10 Au- 
gust 1994 using the Low -Resolution Imaging Spec- 
trometer ( Oke et al. 1995 ) on the 10 m Keck telescope. 
Dithered exposures were taken with total integration 
times of 40 minutes in R and 33 minutes in I. The con- 
ditions were good, with seeing of ~0 / /8 in I and ~ I'll 
in R. These images showed a faint object at the pulsar 
position, but because the field is very crowded, it was 
not clear whether or not the object was the result of 
a blend, and its magnitude was difficult to determine. 

Therefore, the Wide Field Planetary Camera 2 
aboard the Hubble Space Telescope was used to ob- 
serve the PSR B1855+09 held for one orbit each on 
6 August 1997 and 19 February 1998. The first ob- 
servation consisted of three 700 s exposures through 
the F555W filter (mean wa velength A = 5397 A, ef - 
fetive width AA = 1226 A; [Holtzman et al. 19954 , 
the second of five 260 s exposures through F814W 
(A = 7924 A, AA = 1500 A). The field around the 
pulsar was put on a clean spot on the CCD of the 
Planetary Camera (PC). Only the PC images are used 
here. 

For both data sets, the images were registered to 
half-pixel accuracy and resampled using pixels with 
half the original size. Next, cosmic ray hits were iden- 
tified by comparing values at a given position with 
the minimum value occurring among the images at 
that position. A stacked image was formed using all 
unaffected pixels, and this image was resampled to 
the original pixel scale for the further analysis. The 
stacked images are presented in Figure |l|. 

Astro metry was done r elative to the USNO-A2.0 
catalog ( Monet et al. 1998] ), using only those 189 ob- 
jects that were not overexposed in a 10 s I-band Keck 
image and that appeared stellar and unblended. We 
measured their centroids and corrected for instrumen- 
tal distortion using a bi-cubic function determined by 
J. Cohen (1995, private communication). Next, we 



fitted for zero-point position and position angle on 
the sky, keeping the plate scale at the known value. 
Rejecting 28 outliers (residuals > 1"), the inferred 
single-star measurement error is 0'.'33 in each coordi- 
nate. This is somewhat larger than we have found 
for other fields (e.g., Van Kerkwijk fc Kulkarni 1999| ), 
probably because in this crowded held there are resid- 
ual problems with blending in the USNO-A2.0 posi- 
tions. 

The solution was transferred to the F555W and 
F814W PC frames using 50 and 54 transfer stars, 
respectively. For these, the PC positions were cor- 
rected for instrumental distortion using the bi-cubic 



function given by Holtzman et al. (1995b). We fitted 
for the zero-point position only, as the plate scale and 
position angle on the sky are known accurately (see 



Holtzman ct al. 1995b ). Rejecting four and three out- 



liers for F555W and F814W, respectively (residuals 
> 0'.'08), the inferred single-star measurement errors 
for both are 0'.'027 in each coordinate^. 

The main uncertainty in our astrometric solu- 
tion is the extent to which the timing position of 
PSR B1855+09 and the USNO-A2.0 catalog are on 
the same astrometric system. The former is based 
on the DE200 dynamical ephemeris and should be 
close to the International Celestial Reference System 
(ICRS; [Folkner et al. 1994| ). The USNO-A2.0 catalo g 
is tied to the ICRS as well (see [Monet ct al. 1998| ). 
There may be small systematic differences, however, 
as well as effects of the average proper motion of the 
USNO-A2.0 stars between the plate epoch (~ 1954) 
and the time of our Keck observation (1994.6) due to 
Galactic rotation, asymmetric drift, and reflex of the 
solar motion. 

We tried to measure the offset of our solution from 



the ICRS using two Hipparcos (ESA 1997) stars, 
HIP 93083 and HIP 93084, which are on our 10 s I- 
band image. These stars are strongly overexposed, 
but nonetheless we were able to determine good cen- 
troids by fitting all unaffected pixels, with different 
fitting methods agreeing to within 0'.'03. The off- 
sets in both right ascension and declination between 
the positions inferred from our solution and the two 
Hipparcos stars (at epoch 1994.6) were consistent to 
within 0'.'04 (0.02pix), the averages being 0'.'17±0'.'05 
and 0'.'34 ± 0'.'05, respectively. 

For verification, we took 20 AGN with VLBI-ICRS 



5 Leaving scale and position angle free, this reduces to 0'.'014. 
The final position, however, changes by < 0''006. 



3 



43'20" 



05 

I 

o 
o 
o 

0} 



43' 15" 



-i | t^t i i i r 
F555W 



_i L 



_i L 



.'■ r 



• * 



F814W 



« % 



_L 



_L 



57m36s5 



57m36so 

«2000" 18h 



57m36s5 



5711136SO 



Fig. 1— HST images of the PSR B1855+09 field, taken through the F555W (left) and F814W (right) filters. 
The epoch 1998.0 timing position is shown by the tick marks. These are 0'.'24 long, equal to the 95% confidence 
diameter inferred from the uncertainty in the astrometric tie. 



positions within 45° of PSR B1855+09 from the list of 
Ma et al. ( 1998 ) and measured the positions of their 
optical counterparts on the second Digitized Sky Sur- 
vey (epoch ~1991), with astrometry tied to USNO- 
A2.0 (epoch ~1954). We find average offsets consis- 
tent with the above, of 0'.'18 ± 0"04 and 0'.'28 ± 0'.'04 
(the average offsets between the VLBI and the USNO- 
A2.0 positions are 0709 and O'.'ll; the additional dif- 
ference is consistent with what is expected from the 



ave rage proper motion of reference stars at ~ 2 kpc) . 



We conclude that, after correction for the offset found 
using the Hipparcos stars, our astrometry should be 
on the ICRS to 0705. 

Near the position of PSR B1855+09 - aj2ooo = 
18 h 57 m 36!3917, (5j 2 ooo = 09°43'17 7275 at epoch 
1998.0 flKaspi, Taylor, fc Ryba 1994| ) - the PC im- 
ages show one object. It is offset by —0701 ± 0705 
and —0704 ± 0705 in right ascension and declination, 
respectively (see Fig. |l|). Given the density of objects 
of about 1 per three square arcseconds, the probabil- 
ity of a chance coincidence in our 0712 radius 95% 
confidence error circle is ~ 1.5%. Thus, most likely 
we have detected the counterpart of PSR B1855+09. 

Photometry of the object was done following the 
prescription of Holtzman et al. (1995a). We per- 
formed aperture photometry for a range of different 



radii, and used some two dozen brighter stars in the 
frame to determine aperture corrections relative to 
the standard 075 (llpix) radius aperture. For the 
candidate, the best signal-to-noise ratio is for rela- 
tively small apertures, with radii between 1.5 and 
3pix. From these, we infer a count rate for the 075 ra- 
dius aperture of 0.042±0.005 and 0.087±0.009 DN s" 1 
for the F555W and F814W filters, respectively (1 DN 



is approximately 7 detected photons; Holtzman et al 



1995a). These count rates correspond to magnitudes 
of m F555 w = 25.90 ± 0.12 and m F8 i 4W = 24.19 ±0.11 
in the Vega system (using ttifsssw = 22.545 and 
WF814W = 21.639 for a count rate of lDNs -1 in the 
PC, and applying a 0.10 mag aperture correction from 
075 radius to "nominal infinity"; Baggett et al. 1997). 



3. Discussion 

The temperature of the white dwarf can be 
constrained using the measured color, mpsssw — 
"1F814W = 1.71 ± 0.16. For this purpose, we correct 
for reddening using the e stimate Eb~v — 0.5 ± 0.1 
of Kulkarni et al. (1991). This estimate is based 
on CO emission and H I 21-cm absorption mea- 
surements towards PSR B1855+09 and is consistent 
with the range 0.5-0.7 estimated in the general direc- 
tion of PSR B1855+09 from reddening of O-F stars 
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(Neckel, Klare, & Sarcander 1980). The redding cor- 
responds to Af 555 w = 1.6±0.3 and -Ef555W-F81 4W = 

0.65 ± 0.13 flSchlegel, Finkbeiner fc Davis 1998| ), and 
we infer (ttifsssw — wfsi4w)o = 1-06 ± 0.21. This 
intrinsic color corresponds to T c g — 4800 ± 800 K 



if the atmosphere were pure hydrogen (Bergeron 
Saijmon, & Wesemael 1995; using the color trans- 



formations of Holtzman et al. 1995a), but could be 
lower than 4000 K for a mixed helium/hydrogen at- 
mosphere. (Note that the surface gravity expected for 
the companion, \ogg ~ 7.2, is outside the range 7.5- 
8.5 cover ed by the atmospheric models of Bergeron et 
al. [1995], and thus we had to extrapolate; the colors, 
however, are not very sensitive to \ogg.) 

A consistency check on the estimated temperature 
is available using the absolute magnitude. From the 
measured white dwarf mass (Mc = 0.258lg q^Mq; 



Kaspi, Taylor, fc Ryba 1994| ), we infer a ra dius R c = 
O.O21i? (using the models of Driebe et al. |l998| ; the 
result is not sensitive to model details). Combined 
with our estimate of T e g, one infers Mboi — 13.9, and, 
using bolometric corrections tabulated by Bergeron et 
al. ( 1995 ), Mp555\v — 14-3. Correcting for reddening, 
this corresponds to a parallax of 1.0 mas, consistent 
with the timing parallax of 1.1 ± 0.3 mas. 

With T e s and Mc, the cooling age of the white 
dwarf can be estimated using a cooling model. As 
discussed in §[1], the main uncertainty in the mod- 
els is the amount of nuclear burning, which is set 
by the thickness of the hydrogen layer surrounding 
the helium core. Generally, hydrogen layers have 
been assumed to be relatively thin, and little account 
has been taken of a possible dependence on stellar 
mass. For instance, Hansen & Phinney (199S) mod- 
eled helium white dwarfs having hydrogen layers of 
3 x 10" 4 M Q and 1O- 6 M . With their models, one in- 
fers a cooling age for the white dwarf here of ^3 Gyr. 
Using simple evolutionary models for the progenitor 
of the white dwarf, however, much thicker layers are 
found; e.g., for a 1M progenitor tha t ends up as 
a 0.259 Mq white dwarf, Driebe et al. (|l998|) find a 
hydrogen layer of 4.8 x 10 -3 M© at formation. With 
this much thicker layer, nuclear burning is much more 
important, and the cooling age becomes much longer, 
%d = 10±2 Gyr (by which time only ~ 0.7 x 10~ 3 Mq 
of hydrogen is left). 

The white dwarf age inferred from the Driebe et 
al. model is greater than the characteristic age of 
the pulsar, r c = 5 Gyr. If the model were correct, 
t c must be an underestimate of the true age of the 



PSR B1855+09 system. One interpretation is that 
the braking index of this millisecond pulsar is less 
than the canonical value of 3; to obtain tpsR > 8 Gyr 
would require n < 2.25 (Eq. Q). This is perhaps not 
unreasonable, as for most pulsars for which braking 
indices have been measured, values less than 3 have 
been found: n = 2.51 ± 0.01 fo r PSR B0531+21 
flLync, Pritchard, fc Smith 1993j); 2.28 ± 0.02 for 
PSR B0540-69 QBoyd et al. 1995| ); 1.4 ± 0.2 for 



PSR B0833-45 ( |Lyne et al. 1996;); an d 2.837 ± 0.001 



for PSR B1509-58 flKaspi ct al. 1994[ ). All these pul 



sars, however, are young and have strong magnetic 
fields, so their relevance to the discussion here is not 
clear. We note that a variant on the vacuum-dipole 



model ( Melatos 1997 ), which does a reasonable job of 
explaining these braking indices, predicts n = 3 for 
a pulsar like PSR B1855+09 (A. Melatos & J. Hib- 
schman, 1999, private communication). 

The optical counterparts of other pulsar binaries 
may give a clue to where the problem lies. From the 
list compiled by Hansen & Phinney (199S), we find 
that two pulsars, PSRs J0034-0534 and J1713+0747, 
have very cool companions, with T c g < 3500 K and 
T e g = 3400 ± 300 K, respectively. For such temper- 
atures, the cooling ages inferred from the models of 



Driebe et al. (1998) are well over 10 Gyr, even if the 
orbital inclinations were such that the helium white 
dwarfs had close to the maximum mas sQ. While these 
ages are not inconsistent with the estimated pulsar 
ages, they exceed estimates of the age of the Galaxy 
from the white-dwarf luminosity function (for recent 



determinations, 


see 


Leggett, Ruiz, & Bergeron 1998 


Knox, Hawkins, 


& Hambly 199E). This suggests that 



the models may overestimate the cooling ages. 

The above is in contrast to what is the case for 
PSR J1012+5307, where the cooling age estimated 
using the Driebe et al. model is very similar to r c 
(§|l|). The discrepancy might be resolved by the thick- 
ness of the hydrogen layer being a function of the 
orbital separation, perhaps via somewhat different 
mass-loss histories. PSR J1012+5307 has the second- 
shortest orbital period of all systems known (0.6 d), 
much shorter than that of PSR B1855+09 (12.3d). 
PSR J0034— 0534, however, discussed above as an- 



6 The companion of PSR J0034-0534 could also be a CO 
white dwarf, which would make the co oling age much shorter 
( ^ichonberner, Driebe, fc Blocker 1999[ ). The required low or- 
bital inclination has ~ 10% a priori likelihood. The companion 
of PSR J1713 has mass 27 < M c < 0.4 M e and thus shou ld 
be a helium white dwarf (Camilo, Foster, & Wolszczan 1994). 
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other case for which thick hydrogen layers may be 
problematic, has a short orbital period (1.6d). 

The discrepancy might also result from differences 
in white-dwarf mass. For instance, Driebe et al. 



(f999) find that shell flashes only occur in a lim- 
ited range of masses (0.2f-0.30 M©), which includes 
PSR B1855+09 but not PSR J1012+5307. Driebe 
et al. found that burning during the flashes does not 
greatly affect the cooling ages. Schonberner, Driebe, 
& Blocker (1999), however, in duscussing our re- 
sult, have suggested that envelope expansion during 
and following a flash could cause Roche-lobe over- 
flow. The resulting mass transfer to the pulsar (which 
might temporarily become an X-ray source again) 
could lead to a thinner hydrogen layer. If so, then 
with stronger constraints on companion masses, as 
can be derived for PSR JI7I3+0747 in particular, the 
mass range in which shell flashes occur can be con- 
strained observationally. 
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